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ABSTRACT 

We analyse C^®0 (J = 3—2) data from a sample of 99 infrared-bright massive 
young stellar objects (MYSOs) and compact Hii regions that were identified as po¬ 
tential molecular-outflow sources in the Red MSX source (RMS) survey. We extract a 
distance limited (D < 6 kpc) sample shown to be representative of star formation cover¬ 
ing the transition between the source types. At the spatial resolution probed, Larson- 
like relationships are found for these cores, though the alternative explanation, that 
Larson’s relations arise where surface-density-limited samples are considered, is also 
consistent with our data. There are no significant differences found between source 
properties for the MYSOs and Hii regions, suggesting that the core properties are 
established prior to the formation of massive stars, which subsequently have little 
impact at the later evolutionary stages investigated. There is a strong correlation be¬ 
tween dust-continuum and C^®0-gas masses, supporting the interpretation that both 
trace the same material in these IR-bright sources. A clear linear relationship is seen 
between the independently established core masses and luminosities. The position of 
MYSOs and compact Hii regions in the mass-luminosity plane is consistent with the 
luminosity expected a cluster of protostars when using a ^40 percent star-formation 
efficiency and indicates that they are at a similar evolutionary stage, near the end of 
the accretion phase. 
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1 INTRODUCTION 

Massive stars (> SM©) are responsible for some of the 
most energetic phenomena in the Galaxy. They deposit large 
amounts of radiation, kinetic energy and enriched material 
into the interstellar medium (ISM) throughout their forma¬ 
tion, main-sequence lifetimes and when they explode as su¬ 
pernovae. Massive young stellar objects (MYSOs) are the 
precursors to massive stars and are luminous (> 10® L©), 
mid-infrared point source s which have not yet begun to 
ionise their surroundings ll Davies et al. IH). The details 
of their early formation stages are difficult to probe obser- 
vationally, due to their rare, clustered and embedded nature 

* http://rms.leeds.ac.uk/cgi-bin/public/RMS_DATABASE.cgi 
t E-mail:maud@strw.leidenuniv.nl (LTM) 


l|Cesaroni et al ][20 o 3). As a result, the processes of mas¬ 
sive star formation are still relatively uncertain when com- 
pared to the wel l studied low-mass star-formation paradigm 
jshu et al.lll987l ). 

Theoretical modelling suggests that, during their for¬ 
mation, MYSOs with high accretion rates (10~^ M 0 yr~^) 
swell due to the mass influx. Consequently, they are deficient 
in ultraviolet (UV) photons until the y begin contracting to 
a near-Main-Sequence confi guration llHosokawa fc Omukail 
I2OO9I : iHosokawa et al.ll 2010 l) : thus despite being very lumi¬ 
nous, they are not initially ionising their surroundings. As 
the MYSOs contract towards the main sequence, they will 
then start to ionise the surrounding environme nt and form 
expanding Hll regions (|Hoare fc Francol l2007l ). When the 
central stars reach the main sequence they generate copious 
amounts of UV photons, and so rapidly disrupt and destroy 
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the natal cloud. Identifying MYSOs and very young Hll re - 
gions (very compact radio emitters, iLumsden et al.l l2013l l 
provides a sample of sources in which the natal environment 
will be less disrupted and therefore closer to their initial con¬ 
ditions, while simultaneously facilitating the investigation of 
feedback from the massive protostars known to be forming . 
The Red MSX Source (RMS) survey (iLumsden et al.ll2013 l 
is an ideal sample for this since it spans a wide range of lu¬ 
minosity and evolutionary stage. It also allows us to study 
how molecular gas properties vary as a function of both time 
and source luminosity. 

This paper deals with the core properties of a sample of 
primarily northern RMS sources. These observations consti¬ 
tute the only RMS dataset so far in which the molecular-gas 
emission has been mapped, as opposed to obtain ing single¬ 
pointing observations (cf. lUrauhart et al.l|2011^ . The pri¬ 
mary goal of the observations was to study molecular out¬ 
flows using single-dish observations of ^^CO, and 0*^*0 

(J = 3—2) (discussed in a companion paper: Maud et al. 
2015 submitted to MNRAS ). The aim of the 0*^*0 obser¬ 
vations was to study the kinematical behaviour of the gas 
in the molecular core around the RMS sources. This allows 
us to determine outflow properties more reliably, but also 
permits us to study the cores sep arately. Single-dish 
(3—2) is generally optically thin (IZhang fc Gadl2009l l and 
excited in the denser regions when compared with lower- 
excitation transitions of CO (e.g. J = 1—0). The J = 3—2 
emission is much less confused with line-of-sight emission 
than lower-excitation transitions, and a critical density for 
the transition of > lO'^cm”^ (with most emi ssion expected to 
come from significantly higher densities - cf. ICurtis fc Richeij 
l201lh . makes it an ideal tracer of the core structures in star¬ 
forming regions. 

Section [2] describes the source sample and observations 
undertaken, with Section[3]detailing the method used to cal¬ 
culate all source masses and their radii. Section |4] presents 
the basic results and the comparisons with previous data 
from the literature and the RMS survey itself. Section [5] dis¬ 
cusses various relationships with source properties and in 
Section we compare our results with a simple model to 
examine the star formation efficiency and protostellar evo¬ 
lution. A summary is given in Section [7] 


2 SAMPLE AND OBSERVATIONS 
2.1 Sample Selection 

The sources were chosen from all MYSOs and Hll regions in 
the RMS survey that are located within a distance of ~6 kpc, 
have luminosities > 3000 L©, and are observable with the 
JCMT (declinations —25° to -1-65°), with some additional 
right ascension constraints set by the observing dates. In ad¬ 
dition, for the Hll regions, only those sources which appear 
compact in higher resolution mid-infrared images were se¬ 
lected. Finally although all of the sources with L > 10000 L© 
were observed, only a random sample of the less luminous 
ones were included (see below). The original selection was 
made using the pre-2008 RMS catalogue and resulted in 99 
target sources representative of the ~200 MYSOs and com¬ 
pact Hll regions satisfying these criteria in the catalogue 
at that time. Since 2008, the RMS catalogue has evolved 


significantly as a result of subsequent observations. The up¬ 
dates to the RMS database have resulted in 10 of the ini¬ 
tial sources now being assigned a kinematic distance be- 
yond 6 kpc, where lumi nosities are only complete to ^10“* L© 
dMottram et al.lboilal h these are removed from the statisti¬ 
cal sample (although masses are still calculated where possi¬ 
ble). The luminosities for all sources have b een recalculated 
using the m ost up-to-date source distances dUrauhart et al.l 
|2012| . l2014al ') and m ulti-wavelength s pectra l energy distribu¬ 
tion (SED) fits from iMottram et al.l feoilbll . As luminosities 
have changed the sample now extends to ~1000Lq, but is 
not complete below the original limits. Finally, there are 
now a total of ~311 sources in the RMS catalogue that sat¬ 
isfy the original distance-limited criteria described above. 
Of these, we would expect about 50 Hii regions to be too 
extended ( based upon the RMS classifications), so that a 
complete distance-limited sample from the final catalogue 
would number approximately 260 sources. 

The left panel of Figure [1] shows the luminosity distri¬ 
butions of the remaining 89 sources (dark grey) and of ~450 
objects in the RMS database (light grey) now matching the 
original criteria (but where L>1000 Lq and source types are 
MYSO and both compact and extended Hll regions). The 
right panel shows the luminosities as a function of distance in 
the two sets of sources (filled and open circles, respectively). 
Both plots indicate that the range of luminosities and range 
of distances probed by the original observed sample is rep¬ 
resentative of all sources in the RMS database that satisfy 
the selection criteria. It is clear that the observed sample is 
incomplete in terms of relative numbers of sources at lumi¬ 
nosities below 10'* Lq (above which the RMS survey is now 
complete). This was always intended however, since our ini¬ 
tial sampling aimed to have approximately equal numbers 
per logarithmic luminosity bin. Since all the analyses in this 
paper are essentially of ratios of observed quantities (e.g. 
Section where luminosity is concerned), this numerical 
incompleteness is not a significant problem. 

2.2 Observations 

All 99 sources were observed with the James Clerk Maxwell 
Telescope (JCMT) as part of projects M07AU08, M07BU16, 
M08AU19 and M08BU18 during 2007 and 2008. The 15-m 
dish yields a full-width half-maximum (FWHM) beam size 
of 15.3 arcsec at ~329 GHz for the C*®0 (3-2) line. Through¬ 
out the observations, the typical median system tempera¬ 
tures were Tgys ~350—550 K. In a few cases, system tem¬ 
peratures reached as high as 900 K, reflected in a higher 
spectral noise level. The observations were taken with the 
Heterodyne Array Receiver Pro gram (HARP) 16-pixel SSB 
SIS receiver (iBuckle et al.ll200^ . The backend ACSIS corre¬ 
lator (Auto-Correlation Spectral Imaging System) was con- 
hgured with an operational bandwidth of 250 MHz for the 
C*®0 transition. *®CO was simultaneously observed with 
C*®0. The resulting velocity resolution was ~0.06kms“*. 
The C*®0 and *®CO data were re-sampled to the velocity 
resolution of the *^CO data (~0.4kms“*) taken as part of 
the project, in order to match velocity bins for later outflow 
analysis and to improve the signal-to-noise ratio. 

The maps were taken in raster-scan mode with continu¬ 
ous (on-the-fly) sampling and position switching to observe 
a ‘clean’ reference position at the end of each scan row. The 
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Figure 1. Left; Luminosity distribution of RMS sources now conforming to the original sample criteria (light grey bars) over-plotted with 
the luminosity distribution of the distance-limited set of observed sources (dark grey bars). The range of luminosities is well covered. 
Right: Comparison of the RMS distances and luminosities of the 89 sources in the distance-limited sample (filled circles) with the 
MYSOs and Hll regions now conforming to the original sample criteria (open circles) (see text for details). The sample covers the range 
of luminosity in the present population of ~450 RMS sources that now match the original criteria. Note that sourc es at the same distance 
form part of known molecular complexes. Uncertainties are not indicated although are ~ 30 percent for luminosity dMottram et al.l2011hl l 
and an illustrative value of 50 percent can be adopted for the distance. 


maps range in size from 5 square arcminutes (5x5 arcmin- 
utes) for those sources at distances <4 kpc, to 3 square ar¬ 
cminutes for those between 4 and 6 kpc. The pointing was 
checked with reference to a known bright molecular source 
prior to each source observation. Pointing accuracy is within 
~5arcsec, as typically expected from JCMT observations. 
The majority of the baselines on each of the working re¬ 
ceivers were flat, suitable for detecting the weak C^®0 emis¬ 
sion. Some receivers did exhibit sinusoidal modulations and 
were flagged out from the final maps accordingly. Data re¬ 
duction and display were undertaken with a custom pipeline 
which utilised the KAPPA, SMURF, GAIA and SPLAT pack¬ 
ages which are part of the STARLINK software maintained 
by the Joint Astronomy Centre (JAC)0- Linear baselines 
were fitted to the source spectra over emission-free channels 
and subtracted from the data cubes. The final 0*^*0 cubes 
used in all analyses were made with a 7-arcsec spatial pixel 
scale. The data, o riginally on the corre cted antenna tem¬ 
perature scale (r*: lKutner fc Ulichlll98ll ) were converted to 
main-beam brightness temperature Tmb = TX/rimh, where 
Vmb = 0.6 6 as measured by J AC during the commissioning 
of HARP llBuckle et al.ll2009l) and via ongoing planet obser¬ 
vations. Typical spectra noise levels are STmh ~ 0.8 K in a 
0.4kms“^ bin. 


3 MASS AND RADIUS DETERMINATION 

The column density and, hence, mass can be calculated if 
we assume LTE, and a constant Tex and ris for each source. 


^ http://starlink.eao.hawaii.edu/starlink 


as outlined in Appendix A. The calculations rely on the 
accurate determination of J Tmb.is dv over the source. The 
method used here is applicable to any molecular-line tracer 
and source geometry. It is the combined process of integra¬ 
tion over velocity (i.e. creation of a moment zero map) fol¬ 
lowed by an aperture summation over the source area. The 
column density and therefore mass are calculated in each 
pixel via. 


Ar(C^®0) = 5.0 X 10 

X 


12 exp(16.74/Tex) (Tex + 0.93) 
exp(-15.80/Tex) 

f Tmb, 18 71 7 77 dv , 

J [1 - exp(-Ti8)] 


where A(C^®0) is in cm ^ and. 


( 1 ) 


Mgas 


A(C^®0) 


H2 

CiSQ 


fig m(H2)H 


( 2 ) 


where H is the solid angle of a pixel, D is the distance to 
the source, (H 2 /C^® 0 ) is the H 2 to C^®0 abundance ratio, 
where H 2 /'^^CO = lO'^ and is var ied according to 

the relationship 58.8 x T)gc(kpc) -I- 37.1 llWilson fc RoodI 
Il994h . and pg = 1.36 is the total gas mass relative to H 2 . A 
more detailed derivation of the column density and mass are 
given in Appendix|^ (note the constants to change units are 
not included in Equation [2] above). 

The core velocity extent (integration range) is estab¬ 
lished via a direct investigation of the data cubes, channel 
by channel, specifically focusing on a 3-pixel-diameter 
region centred on the source location. The integration 
limits are set when all emission within this region drops 
below (where is the standard deviation from 
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the line-free sections of the re-binned, ~0.4-kms“^ res¬ 
olution spectra extracted at every pixel) while moving 
away from the source Vlsr in the directions of increasing 
and decreasing velocity. For the majority of sources the 
Sctmap contour level is directly traced in order to define 
a polygon aperture selecting which pixels to associate 
with the core, from the moment zero map after veloc¬ 
ity integration (see Figure [H centre, where the dotted 
contour and red contour are Sctmap and the aperture 
respectively). Some sources do not follow this aperture 
definition and are discussed below. All but five sources 
(the YSOs G023.6566-00.1273, G094.3228-00.1671, 

G108.4714-02.8176 and G125.7795-f01.7285 and the Hii 
region G049.5531—00.3302) have strong C^®0 emission 
above the 3 spectral-noise level. These 5 sources are 
not included in further analysis. 

Figure [2] depicts the stages in the process for the 
source G078.1224-f03.6320 (IRAS 20126-f4104, a well stud¬ 
ied MYSO). The left panel shows the average spectrum ex¬ 
tracted in a 3-pixel diameter region centred on the source, 
the integration ranges indicated as dashed lines, and the 
centre panel is the map resulting from the integration in 
velocity. The mass calculated within the dehned aperture 
tracing the 3 u m ap co ntour is ~15OM0, consistent with 
IShepherd et al.l ll2000l . who obtained 104 Mq from their 
G ^®0 ( 1 — 0 ) interferometric observations) when using their 
Galactocentric distance of 8.1 kpc, a heliocentric source dis¬ 
tance of 1.7 kpc, a calculated temperature of 26.1 K and 
without a correction made for ris. As a final consistency 
check, we sum all the data within the aperture into a single 
spectrum and derive a mass from a Gaussian fit to this pro- 
Hle, as shown in the right panel. The result in this case is 
149 ± 6 M 0 , consistent with the previous estimate from ve¬ 
locity integration and aperture summation. Gaussian htting 
of each individual spectrum (at each pixel) in the data cube 
is not used, however, as this would require a priori knowl¬ 
edge of the source size and emission region, which is only 
established after the integration stage. 

The source radii are calculated using the area within the 
defined apertures. An effective circular radius can simply be 
defined as 

0eff = \/ areal'll. 

The deconvolved radius, assuming the JCMT has a Gaussian 
beam of diameter 15.3 arcseconds is 

Sdecor. = 1/2 ([26lefr]" - 15.3")^^" . 

We tested for the influence that extended, low-surface- 
brightness emission might have on these values by clipping 
the data below 30 percent of the peak value rather than 
at the 3 ctmap contour. Although the deconvolved radii de¬ 
crease, the masses also decrease by a similar factor when 
calculated within the same area. Virial-type analyses are 
therefore resilient against the choice of threshold. The mass 
and radius attributed to the cores are therefore consistent 
and refer to the same a rea as in previous similar studies (cf. 
iKauffmann et aI1l2013l i. 

Table [T] lists the sources, positions and important pa¬ 
rameters extracted from the RMS survey online archive. In 
some cases, the types are listed as YSO/Hll where observ¬ 
able characteristics are consistent with both the MYSO and 


Hll-region classification and a d ehnitive type cannot be as¬ 
certained llLumsden et al.ll2013l l. In some sources, multiple, 
close (a few arcsec separation), IR-bright targets have been 
identified (three targets would be listed as A, B, C in the 
online archive, for example). The YSO/Hll type is also used 
for these sources, if at least one MYSO and one Hii region is 
included. Individual MYSOs and Hll regions are inseparable 
at the resolution of the JCMT observations. Furthermore, 
where luminosities for each target have been estimated, the 
total for the source is listed in Table [T] highlighted with an 
asterisk (*). 

Tables [2] and [ 3 ] list the derived source parameters, 
including temperature, optical depths, source radii and 
masses. In 4 of the 94 cores, the ^®CO emission is optically 
thin T 13 <1 and therefore cannot be used reliably to esti¬ 
mate the excitation temperature. For these sources, Tex is 
obtained from the '^^GO data obtained as part of this project 
using T 12 . These sources are flagged with an asterisk (*) in 
Table [2] Although the data taken as part of this work 

is optically thick in all cores, these data are not used to 
establish the temperature of the C^®0 because: the ri 2 =1 
surface is further from the core than that of the which 
more closely matches the region of C^®0 emission; some of 
the spectra also exhibit very strong self absorption; 

and the temperature may be elevated in some sources 

by the influence of outflows. 

Not all of our source maps show an ideal source distri¬ 
bution, such as the test source G078.1224-1-03.6320, Figure 
[ 2 ] Some include diffuse emission or joined targets, thus they 
have a different morphology and their mass apertures do not 
trace the 3crMAP contour directly. In Table [3] each source 
is given a mass flag according to its map morphology and 
aperture dehnition, as detailed in the table caption. These 
flags indicate whether the source masses are reliable esti¬ 
mates and are used in the analysis, or are considered unreli¬ 
able, due to confusion or merging with other sources. Some 
sources are deemed to have unreliable masses where we are 
unable to define a mass aperture following the above out¬ 
lined method, i.e. tracing the 3(7map contour level. Figure[3] 
presents an example source from each flag category, except 
flag types 2 and 6. The masses from sources with 0, 1 and 
2 flags are used in the analysis fully as their mass apertures 
essentially follow the outlined method (see the Table |3] cap¬ 
tion), although flag 4 sources are used with caution where 
indicated throughout the paper. Sources with the other mass 
flags are not used in the analysis as the masses are clear un¬ 
derestimates, the source emission is indistinguishable from 
other cores in the region or emission is not associated with 
the source itself ( see the Table |3] caption in correspondence 
with Figure [3] for the list of mass flags and differences in 
aperture dehnitions). Overall there are 61 (of 89 D < 6 kpc) 
sources flagged as 0, 1 and 2 that are regarded as having 
good mass and luminosity estimates for the cores (70 when 
including flag 4 targets), these all have apertures closely 
following the 3crMAP contour, the remaining sources have 
different aperture dehnitions and are therefore not used in 
further analysis. The integrated images, polygon apertures 
and summed Gaussian spectra of all 99 sources (including 
sources with distance > 6 kpc) are presented in Appendix 
B (available online). The sources with hags 0, 1, 2 and 4 
for which we have derived masses are plotted in all hgures 
henceforth (unless otherwise indicated). 
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Table 1. Source parameters for all objects in the sample, taken from the RMS survey online archive. The asterisk (*) highlights sources 
where multiple IR targets have been identified. Only a small portion of the data is provided here, the full table is available in the 
electronic supplementary information. 


MSX Source Name 

RA. 

(J2000) 

DEC. 

(J2000) 

Type 

Vlsr 

(kms“^) 

Distance 

(kpc) 

Luminosity 

(Lq) 

IRAS source 
(offset) 

Other 

Associations 

GOIO.8411-02.5919 

18:19:12 

-20:47:30 

YSO 

11.4 

1.9 

24000 

18162-2048 (4") 

GGD27 

G012.0260-00.0317 

18:12:01 

-18:31:55 

YSO 

110.6 

11.1 

32000 

18090-1832 (3”) 


G012.9090-00.2607 

18:14:39 

-17:52:02 

YSO 

35.8 

2.4 

32000 

18117-1753 (11^0 

W33A 

G013.6562-00.5997 

18:17:24 

-17:22:14 

YSO 

48.0 

4.1 

14000 

18144-1723 (2") 


G017.6380+00.1566 

18:22:26 

-13:30:12 

YSO 

22.5 

2.2 

100000 

18196-1331 (11^0 


GOIS.3412+01.7681 

18:17:58 

-12:07:24 

YSO 

32.8 

2.9 

22000 

18151-1208 (16”) 


G020.7438-00.0952 

18:29:17 

-10:52:21 

Hii 

59.5 

11.8 

32000 


GRS G020.79-00.06 

G020.7491-00.0898 

18:29:16 

-10:52:01 

Hii 

59.5 

11.8 

37000 


GRS G020.79-00.06 

G020.7617-00.0638=^ 

18:29:12 

-10:50:34 

YSO/Hii 

57.8 

11.8 

62000 


GRS G020.79-00.06 

G023.3891+00.1851 

18:33:14 

-08:23:57 

YSO 

75.4 

4.5 

24000 

18305-0826 (6”) 

GRS G023.64+00.14 


Table 2. Optical depths, main-beam brightness temperatures, excitation temperatures and abundance ratios for all source with strong 
C^®0 emission. Sources with optically thin ^^CO emission are flagged with an asterisk (*) to indicate that ^^CO is used to calculate 
Tex- The full table is available in the electronic supplementary information. 


MSX Source Name 


GOIO.8411-02.5919 
G012.0260-00.0317 
G012.9090-00.2607 
G013.6562-00.5997 
G017.6380-h00.1566 
G018.3412-h01.7681 
G020.7438-00.0952 
G020.7491-00.0898 
G020.7617-00.0638 
G023.3891-h00.1851 


Tnb,13 

+’mb,18 

■^13 


'’"18 


Tex 

Self Abs. 


H2/C'^®0 (10*) 

20.29 ± 0.21 

11.34+ 0.28 

6.07 

+0.35 

-0.33 

0.81 

+0.05 

-0.04 

27.56 12+3 

yes 

7.4 

419.3 

7.78 ± 0.26 

2.42 + 0.34 

2.37 

+0.72 

-0.67 

0.33 

+0.10 

-0.09 

15.21 tlTl 

no 

7.2 

237.0 

13.26 ± 0.75 

7.49 + 0.18 

6.15 

+0.89 

-0.72 

0.83 

+0.12 

-0.10 

20.24 12;|2 

yes 

7.4 

395.8 

7.83 ± 0.30 

3.33 + 0.40 

3.94 

+ 1.00 
-0.86 

0.54 

+0.13 

-0.12 

14.55 

yes 

7.3 

301.7 

14.22 ± 0.20 

8.17+ 0.28 

6.34 

+0.52 

-0.48 

0.85 

+0.07 

-0.06 

91 95 +9.23 

ZL.^O _Q 23 

no 

7.4 

413.4 

20.90 ± 0.24 

7.07 + 0.33 

2.84 

+0.26 

-0.26 

0.38 

+0.04 

-0.03 

29.47 12:22 

no 

7.4 

378.1 

16.70+ 0.24 

5.71 + 0.26 

2.85 

+0.27 

-0.27 

0.39 

+0.03 

-0.04 

24.87 12:22 

no 

7.3 

325.2 

16.70 ± 0.24 

5.71+ 0.26 

2.85 

+0.27 

-0.27 

0.39 

+0.03 

-0.04 

24 87 +9.61 
_Q 52 

no 

7.3 

325.2 

9.58+ 0.23 

2.60 + 0.29 

1.94 

+0.49 

-0.48 

0.26 

+0.07 

-0.06 

18.02 12'22 

no 

7.3 

331.1 

6.32 +0.28 

3.42 + 0.22 

5.67 

+ 1.06 
-0.87 

0.77 

+0.15 
— 0.12 

_o,33 

no 

7.3 

307.6 


Table 3. Measured parameters for all sources. The masses are listed from the polygon aperture method and the Gaussian summation 
fitting where applicable. The full table is available in the electronic supplementary information. 


MSX Source Name 

Type 

Spectral 

Vel. 

Range 

Map Noise 

FWHM 

Centroid 

Aperture 

Gaussian 

Decon. 

Flag 



Noise 






Mass 

Mass 

Radius 




(K) 

(k 

m s“^) 

(K km s“^) 

(kms“^) 

(km s“^) 

(M0) 

(Mq) 

(pe) 


GOIO.8411-02.5919 

YSO 

0.8 

( 9.9 

, 14.6 ) 

1.4 

1.9 

12.2 

465 + 10 

450 + 12 

0.47 

1 

G012.0260-00.0317 

YSO 

1.4 

(110.4 , 

112.1 ) 

1.5 

3.2 

110.9 

269 + 33 

469 + 92 

0.53 

0 

G012.9090-00.2607 

YSO 

0.6 

( 32.6 

, 40.2 ) 

1.7 

4.4 

36.3 

1065 + 79 

1104 + 26 

0.64 

4 

G013.6562-00.5997 

YSO 

1.4 

( 46.4 , 

, 48.9 ) 

1.6 

3.1 

47.8 

108+ 12 

144 + 22 

0.21 

0 

G017.6380+00.1566 

YSO 

0.6 

( 19.9 , 

, 25.9 ) 

1.2 

2.6 

22.3 

616 + 22 

597 + 19 

0.61 

1 

G018.3412+01.7681 

YSO 

0.7 

( 31.1 , 

, 35.8 ) 

1.3 

2.3 

32.8 

606 + 12 

614 + 20 

0.68 

1 

G020.7438-00.0952 

Hii 

1.0 

( 60.0 , 

, 60.8 ) 

0.8 

3.3 

59.1 

594 + 14 

2450 + 150 

1.06 

0 

G020.7491-00.0898 

Hii 

0.8 

( 56.1 , 

, 62.1 ) 

1.8 

3.8 

58.8 

6067 + 147 

6374 + 380 

2.23 

0 

G020.7617-00.0638 

YSO/Hii 

1.0 

( 55.3 , 

, 58.3 ) 

1.3 

3.6 

56.4 

617 + 56 

892 + 101 

0.65 

2 

G023.3891+00.1851 

YSO 

0.8 

( 73.8 , 

, 76.7 ) 

1.0 

2.0 

75.3 

358 + 37 

376 + 30 

0.49 

0 


Mass flags follow the scheme: 

0 — Masses calculated directly from within aperture tracing the 3 ctmap level. 

1 — Faint filamentary structures are not included in mass calculation and are outside the aperture. In extreme cases the aperture 
is more circular. 

2 — Highlights cores with multiple IR-bright sources within JCMT beam (Classically flag type 0 or 1). 

3 — Source mass estimated within a 3-pixel diameter aperture (slightly over 1 beam FWHM) aperture centred on the source due 
to the source being part of a complex filamentary cloud complex. 

4 — Complex/multiple source regions of significant emission. Masses are split where emission peaks are separated by more than 
3 pixels or are circular with a radius set as the shortest distance between the RMS source location and the 3 (Jmap • 

5 — Two or more inseparable continuum cores very close within the aperture. 

6 — Luminosity estimates not from SED fitting. 

7 — Morphology suggesting that gas located at the source position has already been blown away or eroded. 
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Figure 2. Left: 0*^*0 (3—2) line spectrum for G078.1224+03.6320, averaged within a 3-pixel diameter aperture centred on the RMS 
source location and smoothed to 0.4kms“^. The line profile is very close to that of a Gaussian, indicated by the red line. The long 
dashed lines show the velocity range within which the emission is integrated. Centre: Integrated intensity map of the C^®0 emission 
from G078.1224+03.6320 within the —5.7 to —l.lkms“^ velocity range. The dotted black line indicates the 3ctmap contour while the 
solid lines are the contours at 90, 70, 50 and 30 percent of the peak emission value. The white cross indicates the RMS source location. 
Emission is summed within the polygon aperture, shown by the dashed red line tracing the 3(tmap contour level. Right: The summed 
C^®0 spectrum extracted from the data cube within the polygon aperture defined from the integrated map (centre). The resultant 
spectrum has a high signal-to-noise ratio and is Gaussian in shape (red fitted line). The velocity resolution of the summed spectrum 
shown is ~0.1km~^ . 


4 RESULTS 

4.1 Mass and Distance 

Figure H] shows that a wide range of masses are sampled 
at all distances in our sample, indicative of no distance de¬ 
pendent biases. Furthermore, there are no biases towards 
MYSOs and Hll regions independently and the two source 
types span the same mass and distance ranges. In general, 
the sources flagged as 4 also have masses well within the 
range of those identified as having reliable estimates. It is 
likely that the division of mass between sources in multi¬ 
core regions is reasonable in these cases. There are still a 
few flag 4 sources where the masses (>1000 Mq) are likely 
to be overestimates due to the nature of the aperture defini¬ 
tions in these cases, and may include un-associated material 
with the core (see Appendix B). 

All the cores are resolved at distances <6kpc. However 
we cannot resolve sub-structures within the beam. It is likely 
that these cores will form stellar clusters containing lower- 
mass stars as well as the targeted massive protostar (see 
iBontemps et aLll2010l . for an example of a relatively close 
region Cygnus X where substructure is clear). This is sup¬ 
ported by the fact that flag 2 sources are indistinguishable 
from flag 0 and 1 sources in the C^®0 maps, although mul¬ 
tiple infrared sources, with separations less than the JCMT 
beam, are identified in these cores. Clearly we are study¬ 
ing the properties of the natal cores of the associated star 
forming sites here rather than the reservoirs associated with 
individual stars. 


4.2 Comparison with Continuum Masses 

We compared the C^®0 masses with those calculated from 
the 850 Atm integrated flux es from the SCUBA legacy survey 
dpi Fran cesco et aL| 200 ^. BO LOCAM 1.1 mm integrated 
fluxefl (|Cinsburg et al 1 120131) and other 1.2 mm observa- 


^ http://irsa.ipac.caltech.edu/data/BOLOCAM_GPS/ 



Figure 4. Mass versus distance for all sources flagged as 0 and 1 
(circles), 2 (squares) and 4 (triangles). The open symbols repre¬ 
sent MYSOs and the filled ones are Hll regions. There is a wide 
range of masses probed for the different source types and flags lo¬ 
cated at distances <6kpc. Note the error bars in the lower right 
represent a 50 percent uncertainty. 


tions dBeuther et al.l I 2 OO 2 I : iFaiindez et al.l l2004l : iHill et ahl 
[ 200 ^. The continuum fluxes have been used to calculate 
the masses traced by cool dust (Tabled under the optically 
thin assumption via: 


M = 


gS,D^ 

B„{Td) 


(3) 


where Si, is the integrated source flux, g is the gas-to-dust 
ratio = 100, Bi,{Td) is the Planck function for a black- 
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Figure 3. Integrated 0*^*0 (3—2) maps of sources with a range of mass flags. From left to right, top to bottom, are example sources 
indicating the various summation apertures for flags 0, 1, 5, 3, 4, and 7. The source name and corresponding flag is listed above the 
maps. Flag 3 is illustrated by two sources in the G081.7XXX region to exemplify the ~3 pixel diameter summation aperture. For flag 
4, there are 3 sources shown. G012.9090—00.2607 has a circular aperture with a radius set from the source peak to the closest edge of 
the Bitmap noise level. The sources Gill.5320+00.7593 and Gill.5423+00.7776 are a split structure where emission is divided mid way 
between the sources. For G080.9383—00.1268 (flag 7), it is clear that the emission is not peaked on the source location of the Hll region. 
In all maps the dotted black line indicates the Bitmap contour, while the solid contours are the emission at 90, 70, 50 and 30 percent of 
the peak. The white cross indicates the RMS source location and the summation aperture is indicated by the dashed red line. Flag 2 
sources are not shown as they are indistinguishable from type 0 or 1 in the maps, while flag 6 sources are not shown as these can fall 
into any other flag category. 


body at a dust temperature Td and D is the distance to 
the source, k^, is the dust opacity coefficient, calculated via 
Ki^= Kq (vIuq)^, adopting k q = l.Ocm^g "^ at 250 GHz 
( jOssenkopf fc Hennina 1 19941 ) and P — 2 dBeuther et al.l 
l2002h and hence are 1.99 and 1.19cm^g“^ for SCUBA 
(850 pm) and BOLOCAM (1.1 mm), respectively. Assuming 


the gas and dust are in thermal equilibrium, the calculated 
C^®0 gas temperature for each source is used as the effec¬ 
tive dust temperature. This is realistic given the densities of 
such cores (>10'*cm~^. iFontani et ^l2012l ). Furthermore, 
the mean gas temperature for all sources i s i^23K, close to 
typically assumed dust temperatures (e.g. iHill et al]l2005l) 
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Figure 5. Comparison of the C^®0 masses with the 850 pm 
SCUBA masses (circles), the 1.1mm BOLOCAM masses (trian¬ 
gles) and 1.2 mm observations from various other surveys (dia¬ 
monds). There is a dotted line representing equal masses and 
the dot-dashed line indicating the best bisector linear fit of the 
SCUBA and C^^^O masses. These follow each other very closely. 
Open symbols represent reliable C^®0 gas masses (flag 0, 1 and 
2) and the filled symbols are for flag 4 sources (see text). The 
error bars in the lower right represent a 50 percent uncertainty. 


and the kinetic temperatures calcula ted from ammonia ob - 
servations for some of these sources (jUrouhart et al.ll2013h . 
The fluxes listed in the literature are used directly and thus 
the continuum emission regions are unlikely to be exactly 
matched to the emission area of C^®0 and between the dif¬ 
ferent continuum surveys. 

Figure [5] shows that the 850 pm SCUBA, 1.1mm 
BOLOCAM and various 1.2 mm observations correlate 
very well with the C^®0 masses (and each other). The 
BOLOCAM masses plotted here are derived using 80- 
arcsec aperture fluxes which better match the source 
sizes. The error bars represent a 50-percent uncertainty 
in the values. This is a reasonable estimate consider¬ 
ing the calibration accuracy of the fluxes and the poten¬ 
tial variations in dust/gas temperature, dust opacity, in¬ 
tegration aperture selected and the C^®0 abundance ra¬ 
tio (e.g., a factor of ~2.5 difference in mass can be 
caused by a change in as sumed temperature between 10 
and 20K: [Hill et al.ll2005l) . There are a small number of 
outlier sources, G013.6562-00.5997, G025.4118+00.1052, 
G030.8185-b00.2729 and G073.0633-b01.7958 for which the 
dust-continuum masses are noticeably larger than the C^®0 
masses (e.g. > 5 times). Generally, the ~order-of-magnitude 
scatter, may mean the choice of tracer is more important for 
individual objects, if not for the whole sample. However, all 
variations can be explained by a combination of mismatched 
aperture sizes between dust and gas studies, different noise 
levels and the aforementioned temperature and dust opaci¬ 
ties. 

The bisector line of best fit for SCUBA and C^®0 


masses follows the 1:1 line closely in Figured The slopes 
for SCUBA and BOLOCAM fits with the C'^^O masses are 
1.0±0.12 and 1.14±0.12 respectively. The SCUBA observa¬ 
tions appear to trace the same material for the majority 
of these cores. For the cores with largest gas masses the 
BOLOCAM masses are typically greater. However, this is 
likely to be due to the use of a constant aperture size for 
BOLOCAM fluxes and could include more faint, extended 
emission of the more massive regions in comparison to the 
polygon apertures chosen here. This is consistent with the 
association of C'^^O with denser gas as discussed below in 
Section m 


Table [S] presents the results of a formal correlation anal¬ 
ysis for the gas and dust masses for different source types 
and flags. Overall, the method used to calculate the C^®0 
masses produces values that are directly proportional to 
continuum-based mass estimates. At the resolution of the 
JCMT, the C^®0 and SCUBA 850-/im emission effectively 
traces the same regions in the majority of the sources, and 
only in a few cases (where sensitivities are notably different 
between the sets of data) do the emission regions vary signif¬ 
icantly. The BOLOCAM masses calculated from 80-arcsec 
aperture fluxes also provide reasonable matches to the C'^^O 
gas masses but fixed apertures are not ideal for such sources. 
Continuum and C'^^O line emission clearly trace the same 
material for these IR-bright sources. 


In previous studies, depletion of C'^^O is a known cause 
of reduced gas column de nsity and, hence, reduced masses, 
for a range of cores (e. g. ICaselli et al.l 1999yFontaniet_^lJ 
|2012| : lYildiz et al.ll201^ . Table 6 from Bergin et al. ( 19951) 
shows that the amount of CO in the gas phase varies by a 
factor of ~2 between dust temperatures of 20 to 24 K, which 
closely match es the calculated temperatur es of some sources 
in this work. iLopez-Senulcre et al.l (120101) also discuss how 
their gas masses are on average lower than those calculated 
from dust emission although they suggest an incorrect abun¬ 
dance ratio can account for this. Clearly adopting a different 
H 2 /^^C 0 ratio other than 10"^, by a factor of 2, i.e. setting 
H 2 /^^C 0 ~ 2x10“^ will already alleviate any discrepancy in 
dust and gas masses. Factor of 2 offsets have be en found 
in massive, pre-stellar cores (iFontani et al.l l2006l) and at¬ 
tributed to depletion. Here the reasonable correspondence 
of the gas and dust masses indicates that depletion of CO 
is not significant in the majority of these IR-bright RMS 
sources, especially when contrasted with the IR-dark cores 
in the aforementioned studies. Depletion could cause some 
of the scatter we see, however the already discussed vari¬ 
ations in apertures size, dust opacity and temperature for 
example, can also explain this. 


Finally, we stress that the masses calculated from 
the C^®0 emission are homogeneous in their determina¬ 
tion and definitively associated with only the targeted 
cores, given that they have a single velocity component. 
iKauffmann et aP (l2013l) for example, include only cores 
where single Gaussian components are found to avoid the ar¬ 
bitrary division of continuum masses between velocity com¬ 
ponents. These C'^^O observations therefore provide an ideal 
way to study the core mass and its relationship to other ob¬ 
servables. 
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Table 4 . Continuum masses for a number of sources in the sample Table 5 . Table of Spearman rank correlation values comparing 

from BOLOCAM, SCUBA and other 1.2 mm observation. Masses dust and gas masses. The significance of a given p value de pends 

are reasonably consistent between continuum observations. De- on the sample size (see Table A2.5, of lWall Jenkinj|20Q3l . from 

pendent upon dust opacity, temperature and observable uncer- which we have also adopted the quoted P-values). P-values of 

tainties a generous error estimate of ^^50 percent can be assumed 0.05, 0.002 and <0.001 represent the ~2, 3 and >3cr confidence 

(as Figure [Sj. levels. 


MSX Source Name 

Type 

BOLOGAM 

(M©) 

SCUBA 

(M©) 

1.2 mm 

(Mq) 

GOIO.8411-02.5919 

YSO 


139 

210 

G012.0260-00.0317 

YSO 

2188 

3109 

1853 

G012.9090-00.2607 

YSO 

838 

1167 

1528 

G013.6562-00.5997 

YSO 

1142 

1385 

1259 

G017.6380+00.1566 

YSO 

358 

374 


G018.3412+01.7681 

YSO 


224 

307 

G020.7491-00.0898 

Hii 

4016 



G020.7617-00.0638 

YSO/Hii 

1996 



G023.3891+00.1851 

YSO 

1087 



G023.6566-00.1273 

YSO 

584 



G023.7097+00.1701 

Hii 

2569 


2966 

G025.4118+00.1052 

YSO 

886 

1329 

1446 

G028.2007-00.0494 

Hii 

3099 

4783 

3693 

G028.2875-00.3639 

Hii 

8234 

9053 

11894 

G028.3046-00.3871 

YSO 

2537 


1377 

G030.1981-00.1691 

YSO 


372 


G030.6877-00.0729 

Hii 



1956 

G030.7206-00.0826 

Hii 

5666 


2308 

G030.8185+00.2729 

YSO 

432 

440 

415 

G033.3891+00.1989 

YSO 

357 



G037.5536+00.2008 

YSO 

4714 

3552 

3113 

G043.9956-00.0111 

YSO 

408 

597 


G045.0711+00.1325 

Hii 

904 

268 


G048.9897-00.2992 

YSO/Hii 

1621 

1830 


G050.2213-00.6063 

YSO 


397 


G053.9584+00.0317 

Hii 

332 



G073.0633+01.7958 

YSO 


70 


G075.7666+00.3424 

YSO 

138 

88 


G077.9550+00.0058 

Hii 

13 



G077.9637-00.0075 

Hii 

9 



G078.1224+03.6320 

YSO 


90 

138 

G078.8867+00.7087 

YSO 

391 

611 


G079.1272+02.2782 

YSO 


24 

47 

G079.8749+01.1821 

Hii 

54 



G080.8624+00.3827 

YSO 

51 

82 


G080.8645+00.4197 

Hii 

89 

137 


G080.9383-00.1268 

Hii 

31 

18 


G081.7133+00.5589 

Hii 


367 


G081.7220+00.5699 

Hii 

434 

312 


G081.7522+00.5906 

YSO 

260 

272 


G085.4102+00.0032 

YSO/Hii 

1448 

2492 


G094.6028-01.7966 

YSO 


845 


G103.8744+01.8558 

YSO 


91 

117 

G105.5072+00.2294 

YSO 

480 



G105.6270+00.3388 

Hii 

661 

631 


G109.0775-00.3524 

YSO 

404 

179 


G109.0974-00.3458 

Hii 

220 

161 

275 

G109.8715+02.1156 

YSO 

88 

112 


GllO.0931-00.0641 

YSO 

650 

765 

1013 

GllO.1082+00.0473 

Hii 

1008 

1027 


Gill.2348-01.2385 

YSO 

301 

235 

391 

Gill.2552-00.7702 

YSO 

316 

371 

477 

Gill.5234+00.8004 

YSO 

135 

109 


Gill.5320+00.7593 

YSO 

1033 

1119 


Gill.5423+00.7776 

Hii 

1215 

1167 


Gill.5671+00.7517 

YSO 

657 

952 


Gill.5851+00.7976 

YSO 


12 


G133.6945+01.2166 

YSO/Hii 

383 

613 


G133.7150+01.2155 

YSO 


794 


G133.9476+01.0648 

Hii 

572 

1148 


G134.2792+00.8561 

YSO 

40 

64 


G136.3833+02.2666 

YSO 

109 

164 


G138.2957+01.5552 

YSO 

206 

250 


G139.9091+00.1969 

YSO/Hii 


163 


G141.9996+01.8202 

YSO 


66 


G192.5843-00.0417 

Hii 

293 

291 

220 

G192.6005-00.0479 

YSO 

136 

130 

89 

G196.4542-01.6777 

YSO 


279 


G203.3166+02.0564 

YSO 

52 

42 


G207.2654-01.8080 

YSO/Hii 


172 



Correlation 

Size 

P 

P-value 

SCUBA 

All (flag 0+1) 

39 

0.60 

< 0.001 

All (flag 0, 1, 2+4) 

46 

0.60 

< 0.001 

YSO (flag 0, 1, 2+4) 

32 

0.53 

0.001 

Hii (flag 0, 1, 2+4) 

14 

0.64 

0.01 


BOLOCAM 


All (flag 0+1) 

33 

0.50 

0.002 

All (flag 0, 1, 2+4) 

39 

0.58 

<0.001 

YSO (flag 0, 1, 2+4) 

23 

0.22 

0.3 

Hii (flag 0, 1, 2+4) 

16 

0.86 

<0.001 



Figure 6. Comparison of the C^®0 linewidths with NH3 (1,1) 
(filled circles), NH3 (2 ,2) (red triangles), NH3 (3,3) (squares) from 
lUrguhart et alJ ll201ll') .The dotted line is the 1:1 line and is be¬ 
tween the NH3 (2,2) and NH3 (3,3)transition. An average error 
of ~ 5 percent from the fitting procedure is indicated. 


4.3 Comparison with Other Linewidths 

Previous observations of molecular gas, both in our own and 
other galaxies, have led to differences in derived properties 
depending on the species involved, and in particular when 
CO is compa red to tracers of dense molecu lar gas (e.g. the 
discussion in iKrumholz fc Thomps^l2007l) . We have com¬ 
pared our 0*^*0 linewidths ( measured using integrated spec¬ 
tra extracted from within the mass aperture regions, see 
Figure [ 2 ] and Appendix B) wi th those from NH 3 da ta from 
the single-dish observations of lUrauhart et aP ll201ll) in Fig¬ 
ure |6l The increase in linewidth with excitation energy in 
NH 3 is already well reported in similar types of sources (e.g. 
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iLongmore et aL I I2OO7I : IWienen et al.ll2012| ~l. and is generally 
interpreted as due to the hotter gas being nearer the central 
exciting source, and wit h a smaller beam fillin g factor in the 
NH3 observations (e.g. lUrauhart et al.ll201lh . The widths 
are primarily driven by motions within the gas rather than 
temperature broadening, which would give far smaller val¬ 
ues. Our data show that the C^®0 J — 3 — 2 linewidths are 
consistent, in magnitude, with values between the NH3(2,2) 
and NH3(3,3) linewidths. The optically thin 0*^*0 3-2 tran¬ 
sition is tracing similar density (~ 10"^ cm“®) core material 
as the ammonia line, but should trace cooler gas than both 
NH3(2,2) and NH3(3,3), since the rotational energy levels (in 
temperature units) of the CO J = 3 — 2 line (~33 K) and 
NH3(1,1) (~23K) are actually better matched (cf. NH3(2,2) 
~65K and NH3(3,3) ~125K). This suggests that the CO 
ladder is actually thermalised to a higher J level, and that 
the kinematics of the J = 3 — 2 transition in our sources are 
actually representative of that warmer gas as well. We also 
compared the line centre velocities of the NH3 and C^®0, 
and found they agreed well, with no deviation larger than 
the typical spread in the various NH3 linewidths shown in 
Figure [6l Again, this demonstrates that both these tracers 
are sampling similar gas properties and volumes. 

We also compared the J = 3 — 2 C^®0 and ^®CO 
linewidths, with the latter larger by at least 10 percent. The 
largest differences are seen for those objects where the line 
opacity for ^®CO is largest, so this is generally caused by 
the opacity of the CO lines rather than the fact that the 
^®CO line also traces outflow material. A similar trend to¬ 
wards increasing linewidth with increasing line opacity is 
also seen when comparing the C^®0 J = 3 — 2 data with 
our previous lower J transition data used in deriving kine¬ 
matic distances fe.g. lUrauhart et al.l 1200^ 1. Indeed other 
similar surveys have commented on the di fference in CO 
linewidth for the low lying isotopolog ues le.g. lAo et al.ll2003 : 
IDu fc Yanel200^ : IWieneii et al.l20ld l. The key message here 
is that different tracers may be more suitable in different cir¬ 
cumstances. For example, the J = 3 — 2 C^®0 studied here 
is clearly a good tracer of kinematics in reasonably dense 
clumps/cores, whereas J = 1 — 0 ^®CO is more suitable 
for probing the diffuse cool gas that delineates molecular 
clouds as a whole. This also indicates why the J = 1 — 0 
observations are a good tracer of gas in other galaxies, since 
the bulk of the mass will tend to be in the more diffuse 
molecular clouds. Notably the same may not be true in the 
dense environment s found in extreme star forming galaxies 
dHarris et al. IH). 

The mean FWHM values for MYSOs and Hll regions 
are ~2.6±0.1 and ~3.1±0.2kms“^ respectively, where the 
uncertainties are the standard errors (including all sources 
where D<6kpc and C^®0 is detected). The slightly larger 
FWHM for Hll regions might be interpreted as an evolu¬ 
tionary trend, where linewidths increase as the source has a 
greater impact upon its sur roundings. How e ver, s uch an in¬ 
terpretation is incorrect. As lUrauhart et ^ (l201ll ') note, any 
such trend is artificial and caused by the luminosity-FWHM 
relationship due to the different luminosit y functions for 
MYSO s and fully developed Hii regions Isee lMottram et al.l 
l2011bl '). Furthermore, the Hll regions here, as explicitly 
noted in Section [2 were selected to be compact and there¬ 
fore should be at a similar evolutionary stage to the MYSOs. 
Figure [7] shows the relationship between the source luminos- 



Figure 7. Plot of luminosity versus 0*^*0 FWHM for all sources 
flagged as 0 and 1 (circles), 2 (squares) and 4 (triangles). The 
open symbols represent MYSOs and the filled ones are Hll re¬ 
gions respectively. A linear log-log relationship is seen between the 
parameters. The uncertainty from FWHM fitting is ~ 5 percent 
where as that for the lumin osity is a representative 30 percent 
value jMottram et alj|2011bh . 


ity and FWHM. Evidently, the Hll regions have a greater 
proportion of sources at luminosities > 10'* Lq (and corre¬ 
spondingly larger linewidths) which explains the slight offset 
of mean FWHM values reported. The Spearman rank cor¬ 
relation coefficient is 0.45 at <0.01 significance level, inter¬ 
preted as a strong correlation for the 61 sources (with mass 
flags 0, 1 and 2). The relationship does hint at the most 
luminous sources providing more feedback and turbulence 
to the cores, possibly by driving more powerful outflows. 
This is examined further in Maud et al 2015 (submitted to 
MNRAS). 

It is worth noting that some sources (30, 32-including 
those where D>6kpc) show evidence of regular velocity gra¬ 
dients spatially across their cores. However, only a few of 
these appear to be aligned with, or perpendicular to the 
outflow direction. We stress that, given the resolution of the 
observations, it is difficult to ascertain the cause of the gra¬ 
dients in most cases, and whether the outflows are a major 
contributor. 

In order to examine whether the CO lines profiles reveal 
the presence of infalling motions we derive d the line asym¬ 
metry parameter using the definition of iMardones et al.l 

(fi^ . 

vithick) — vithin) 

ov = -T— . , . . - 

Av(tnm) 

where v(thick) and v(thin) are the velocities of the peak pix¬ 
els measured from the spectra and Av{thin) is the FWHM 
of the thin tracer. Here C^®0 is the optically thin tracer (see 
Table [21), and the line is the optically thick one (also 
see Maud et al. 2015 submitted to MNRAS). The optically 
thin FWHM values of C^®0 are those given in TableO There 
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is a spread of asymmetry values, ra nging: from —1 to +1.5 
(where by the lMardones et al.lll997l criterion |5ti| > 0.25 in¬ 
dicates asymmetry). The results show that 23 sources have 
a red asymmetry, and 15 a blue asymmetry. At face value 
this can be taken as a lack of evidenc e for onRoin R infal l 
in our sample. However, we note that [Fuller et al.l (l2005l l 
also find a similar result in the profiles of the transitions 
they study which have similar excitation temperature and 
critical density, even though the same sources show strong 
blue asymmetry in, for example, HCO’^ J = 1 — 0. Our re¬ 
sults appear to add evidence to their comment that it is the 
dense gas tracers that are best tracers of large scale infall. 
CO appears to be a poor tracer of infall for these cores, 
while HCO "*" (4—3) is usually much better suited to these 
size scales (|Klaassen fc WilsorJl2nn7t l. 


5 ANALYSIS 

5.1 Correlations within the Observables 


We examined possible correlations amongst the observable 
parameters present in our sample, and within the wider in¬ 
formation held as part of the RMS database. In addition we 
also considered combinations of these parameters. In par¬ 
ticular we used the derived values from the C ^®0 observa¬ 
tions for radius, Tex, ria, mass, line asymmetry and FWHM 
as well as combinations of these such as gas surface den¬ 
sity, as discussed below. We took information on the lumi¬ 
nosity, infrare d colours and galacto centric radius, Rgc from 
the database llLumsden etah 2013ll. as w e ll as a selection of 
the NHa properties from Urguhart et al.l (1201 il l. We used a 
Kendall correlation method for all of these comparisons. The 
Kendall method is generall y held to be better in the pres¬ 
ence of errors in the data ( Wall fc ,Tenkinsll2003l l. however 
typically both the Kendall and Spearman rank correlations 
lead to the same inferences. Here we choose the Kendall cor¬ 
relation since some parameters have ties in the data, where 
both pairs have the same ordinal. 

The net result is that there are relatively few strong 
correlations, and many of those that are present are driven 
by the one obvious dominant relation, that between mass 
and radius. A partial correlation analysis conhrms this. The 
only other strongly signihcant parameters are luminosity 
and FWHM. Table [S] shows the resultant correlations that 
have signihcance values < 0 . 002 . 

The correlations that are not dependent on mass, lu¬ 
minosity, radius or FWHM can be summarised briefly as 
follows: 


• The NH a kin etic temperatures, Tkin, from 
I Urguhart et al.l ((loTJ) agree reasonably with the Tex 
values we derive here (See Appendix A). 

• The kinetic temperatures are correlated with both in¬ 
frared colour measures we use, Fws/Fk and Fs/Fk, with 
Tex showing a slightly weaker correlation than Tkin (formally 
not significant at our threshold level for Fs/Fk)- 

• ri 3 is anticorrelated with Rgc. 

The hrst of these essentially shows that the CO and NH 3 
are tracing similar material, as we argued previously in Sec¬ 
tion in The second is curious since we might expect the 
redder objects to be more embedded and have lower kinetic 


Table 6 . Table of Kendall rank correlation values for various 
derived observables from both this paper and the wider RMS 
database. R is the deconvolved source radius, M the mass, L the 
luminosity, Av the FWHM, p the density, S the gas surface den¬ 
sity, T 13 the optical depth of *^^00, Rgc the Galactocentric radius, 
the free-fall time, Tex (CO) the excitation temperature of the 
^^CO, TkinfNH s ) the kinetic temperatures derived from NH 3 by 
lUrquhart et alJ ll201lh . Fg/Fj^ the ratio of 8 /rm MSX flux with 
K band flux as derived for iLumsden et al.l ll2013h and Tws /Ffc 
the ratio of the WISE band 3 flux (approximately 10pm) with K 
band flux. The — indicates data where there are ties so no well 
defined significance exists. The significance is the probability that 
the null hypothesis, i.e. there is no correlation, is correct. 


Correlation 

Sample Size 

T 

Significance 

R-M 

58 

0.78 

< 0.0001 

R — Ai; 

58 

0.29 

0.002 

R-L 

58 

0.47 

< 0.0001 

M — Av 

58 

0.39 

< 0.0001 

M - L 

58 

0.51 

< 0.0001 

Av — p 

58 

0.35 

0.0001 

An - S 

58 

0.35 

0.0001 

ns - p 

58 

0.35 

0.0001 

ri3 — Rgc 

58 

-0.37 

< 0.0001 

ns - S 

58 

0.35 

0.0001 

tfj — S 

58 

-0.42 

< 0.0001 

Tex (CO) - Tki„(NH3) 

41 

0.52 

- 

Tkin(NH3) - Fs/Fk 

37 

0.42 

- 

Tki„(NH3) - Fws/Fk 

31 

0.41 

- 

SSFR = S/tfj 

58 

0.49 

< 0.0001 


temperatures rather than higher. It may indicate that these 
warmer sources simply have a more centrally concentrated 
mass distribution. This is not something that we can test 
with the spatial resolution of the current data. Finally the 
anticorrelation between T 13 and Rgc may simply be a re¬ 
flection of the obse rved anti-correlation b etween luminosity 
and Rgc found by iLumsden et al.l ll2013l) in the full RMS 
sample (since T 13 is correlated with density, which is weakly 
correlated at the P= 0.025 significance level with luminos¬ 
ity), even though no significant correlation between these 
variables is found in the much smaller sample here. 

Mass, luminosity, size of clump and FWHM are all pos¬ 
itively correlated with each other (the luminosity-FWHM 
correlation has a significance level of only 0.003 however). 
These four colle ctively are wha t is expected from the scaling 
relationships ofE arsoni lll98ll) . Indeed, appropriate projec¬ 
tions of the mass-luminosity-radius plane, for example, re¬ 
sult in most data points being strongly clustered (i.e., there 
is a mass-radius-luminosity fundamental plane), though the 
luminosity component of this relationship is small, and the 
same is true for mass, FWHM and radius. The first of these 
can be understood, for example, as Larson’s third relation 
(mass-radius) but modified for the case of cores where a not 
insignificant fraction of the mass is now locked into stars 
(and hence not traced by the gas mass: cf. Section | 6 . 1.211 . 
The weak luminosity dependence of this mass-radius rela¬ 
tion effectively compensates by directly tracing the mass 
already locked into stars. The strong correlation of mass 
and luminosity reflects the fact that these sources are of a 
similar evolutionary s t ate, o therwise more scatter would be 
evident. iDavies et al. liEll) shows that Hll regions can be 
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of a similar age to massive protostars if the central source is 
itself more massive (i.e. more massive objects evolve more 
quickly towards the zero age main sequence). This conclu¬ 
sion is justified further in Section It is notable that the 
FWHM and radius (Larson’s hrst relation) gives the weak¬ 
est correlation, as shown in Table [5] We will discuss this 
and the other Larson relations in detail in the following 
sections, as we l l as co mp aring our results with those of 
lUrauhart et al.l ll2014bl ') and iHever et al.l (l2009l 'l . iHever et al.l 
1 2 OO 9 II used J = 1 — 0 data from the Galactic Ring 

Survey (Ijackson et al.llioOfill to re-examine the underlying 
physical principles that lead to Larson’s scaling relations 
for molecular-cloud-sized structures and, hence, sampled 
much larger reg i ons bu t with much lower surface densities. 
lUrouhart et all (|2ni4bh analysed ATLASGAL sources (in¬ 
cluding those with RMS counterparts) and hence studied 
the continuum dust emission of similar regions to us. 

Some of the properties that show no correlation are also 
worth mentioning. The line asymmetry parameter shows 
no correlation with any other properties, including the line 
opacity itself for example. There are no correlations between 
infrared colour and properties such as mass, or line asym¬ 
metry, where both red and blue asymmetric data show the 
same average colour, or between most of the GO properties 
(including opacity and column density) and infrared colour. 
The beamsize of our JCMT data are not dissimilar to many 
of the mid- and far-infrared data we use the infrared colours 
from. The natural explanation for the lack of correlation 
therefore is that there is substructure within the beam for 
both sets of data. In particular the CO data presented here 
traces cooler gas, whereas the infrared colours we have are 
predominantly a combination of hot and warm dust and ex¬ 
tinction. The spatial scales of these components should differ 
considerably but we are unable to probe such detail. This 
also tallies with the discussion above regarding the actual 
correlation seen between Tex or Tkin and colour. 


5.2 Mass and Radius 


Previous observations of molecular clouds/cores suggest a 
power-law scaling between mass and length scale (size, 
radius) of the form m oc r'’', where 7^2 (e.g. 


LarsonlllQSlUElrnegreen fc FalgaroMll996l : iKirk et ^l20ld : 
Kauffmann et alT 2013h . Their observations include differ¬ 
ent regions, associations (low- and high-mass stars) and ob¬ 
servations using both molecular line tracers and continuum 
emission. Figure [ 8 ] shows the G'^^O masses plotted against 
deconvolved radius for our sample. The linear trend when 
including the flag 0 , 1 and 2 sources follows the power law 
m oc where the uncertainty only accounts for the 

spread in the raw mass and radius values. Note, exchanging 
the deconvolved radii for the effective radii has a minimal 
effect and slightly steepens the slope to m oc r -2 i±o i 

The ATLA SGAL clumps a s sociate d with RMS sources 
investigated bv lUrauhart et ^ (l2014lJ) span a larger range 
of radii and mass, and follow a slightly shallower slope 
(7 ~ 1.75±0.04) overall. There is a turnover in their data at 
higher masses, such that if we restrict the range to that of 
our sample a slightly steeper slope ~1.9 would be recovered. 
However, the discrepancy in slopes can be attributed to the 
different methods in which mass and radius are calculated 
(we see that the slope changes when using effective radius). 


lUrauhart et al.l (l2014bll find that the large scale clump prop¬ 
erties follow the same trend from embedded maser sources 
to extended Hll regions (i.e. less to more evolved). They sug¬ 
gest the clump properties must therefore be set prior to the 
onset of star formation and that the subsequent evolution 
of massive protostars, and their feedback, does not effect 
the clumps overall. Note tha t the internal structure s may 
still evolve (cf. the study of iKauffmann et al.l I 2 OI 0 I which 
examined mass as a function of radius within single clumps, 
as opposed to the inter-clump com parisons given here). We 
basically hud the same relations as Urqul 2 jart_et aljj 2014bli 
from the gas as opposed to dust. iHever et al.l 1 2 OO 9 I 1 find 
7 ~ 2.17 ± 0.08, again in reasonable agreement with the 
slope found here. 

Where the Heyer data differs from the smaller 
clump/core sized regions studied here is in the offset of 
this relationship, which lies at lower masses in the Heyer 
et al sample. Fundame ntally, both our sources and those of 
lUrauhart et al.l ll2014b[l are significantl y more massive tha n 
the similar size regions in the sample of IHever et al.l (I2OO9II . 
It is worth reiterating that Heyer et al. used J = 1—0 
data. The comparison in Section 113] for our objects sug¬ 
gests such lines should be broader (due to opacity) than the 
J = 3—2 C^®0. Therefore, if anything, the true offset from 
a fair comparison with the Heyer et al. data may be even 
lower than suggested by the raw data. This is contrary to 
the spirit of the relation that Larson initially proposed. 


5.3 FWHM-Mass and FWHM-Radius 
relationships 

Figure m shows the log-log relationship between the G^®0 
FWHM and deconvolved source radius. This comparison is 
analogous to the type- 2 , single tracer m ulti-core relation¬ 
ships presented in iGoodman et all (Il998l ~) where the FWHM 
is observed to decrease for smaller cores. The meaningful re¬ 
sult of such a relationship is to examine whether the cores 
are virialised, where FWHMoc ® is expected. The bisector 
ht best represents the data given the scat ter (compared wit h 
ordinary least squares fitting, OLS, see llsobe et al.|[l99Cll L 
resulting in FWHMocr° ®^°'®, roughly consistent with viri¬ 
alised cores. However, as already noted (Table O , the sig- 
nihcance of the correlation seen here is very weak (it is no¬ 
tably stronger for the data from IHever et aHl2009ll . As with 
the results of the previous section, a noticeable offset exists 
in the relationship see n here and that of molecular cloud 
scales from the data of IHever et al.l (l2009l L The cores here 
exhibit much larger FWHM values. One possible cause for 
the weaker correlation and larger values is that on these 
scales feedback has a more signihcant impact compared to 
gra vitational moti ons. 

iLarsonl lll98ll ') noted a clear relationship between mass 
and velocity dispersion (linewidth) over many orders of mag¬ 
nitude in mass. Recent studies have conhrmed this, though 
the tight relationship Larson found is less well reproduced. 
As the core masses are strongly related to the source radii, 
which in turn are weakly correlated to the linewidth, there is 
expected to be a link between mass and FWHMs for these 
sources. Figure [10] shows the correlation present between 
the mass and FWHM linewidths. A bisector best ht to the 
data indicates FWHM oc although the OLS ht 

is noticeably shallower, FWHM oc and closer to 
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Figure 8. Core mass versus the deconvolved effective radius as¬ 
suming spherical cores. Open and filled symbols are MYSOs and 
Hii regions respectively. The sources with 0 and 1 flags are rep¬ 
resented by circles, flag 2 by squares and flag 4 sources by tri¬ 
angles. The general linear fit trend follows a power law (7 = 
2.0 di 0.1) indicated by the dot-dash ed line. The dotted line is 
that from iKauffmann Pillail ll2010l) where massive star form¬ 
ing regions exceed a mass of 580i?p{^. The upper and lower light 
grey boxed regions rep resent the high and low m ass star form¬ 
ing regions presented in iKauffmann Pillail l|201Cl) after a reduc¬ 
tion in mass by 1.5 to correct for the difference in dust opacity 
used. Note the larger, more massive ATLASGAL sources from 
lUrguhart et al.l ll2014bh are predominantly located in the HMSF 
region where M^IOOOMq and R>0.5pc. A representative uncer¬ 
tainty of 50 percent in mass and an 10 percent illustrative uncer¬ 
tainty for the radius is indicated to the lower right. 


that found bv iLarso 3 ||1981D who fit by eye. The slopes are 
consistent given the uncertaintie s as a result of the scat¬ 
ter in the data. The lHever et al.l (l2009f) data give a similar 
slope, oc but again the relationship is offset, with 

smaller linewidths for molecular-cloud-sized regions at the 
same masses as in our sample. 

5.4 Virial Mass and Gas Surface Density 

Star formation has already occurred in these cores as they 
harbour at least one IR-bright protostar. However, the virial 
masses can also be used as an additional test to investigate 
the impact of feedback from these sources. Figure [ 11 ] shows 
the core masses against t he calculated virial masses 

following iMacLaren et al.l lll988h : 

M„i,(M 0 ) = 126i?(pc) FWHM^{kms~^), (4) 

where R is the deconvolved radius, we use the FWHM 
of the 0*^*0 emission and assume spherical cores with a 
p oc r~^ density distribution and no magnetic support. 
Changes in geometry, or the density law (Figure [ 8 ] would 
imply a p oc r~^ for a spherical geometry for example), 


Figure 9. Deconvolved source radius versus C^®0 FWHM for 
MYSOs (open) and Hll (filled) regions flagged as 0 and 1 (circles), 
2 (squares) and 4 (triangles). The bisector fit best representing 
the data is shown as the dot-dashed line, and obeys FWHM oc 
j,0.8±0.5 average uncertainty of 5 percent in FWHM and an 
illustrative 10 percent uncertainty in radius is indicated to the 
lower right. 



Figure 10. Core mass versus the FWHM linewidth for MYSOs 
and Hii regions (open and filled symbols respectively). Flag 0 
and 1, 2 and 4 are represented by circles, squares and triangles. 
The plotted dotted line is the OLS best fit with a s i milar slope 
(FWHM oc m,'^ l8±0.04^ by_eye fit from iLarsonI lll98j) . The 

bisector fit better represents the data given the scatter and has 
a steeper slope, FWHM oc . The average uncertainty 

of 5 percent in FWHM and representative 50 percent uncertainty 
for the mass is shown to the bottom right. 
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generally increases the virial mass by up to 50 percent (cf. 
iKaufFmann et ^ .I 2 OI 3 II . Furthermore, we note that the virial 
masses would have artificially been elevated if optically thick 
tracers such as and were used for this analysis 

given their larger linewidths. Using the (3—2) data 

would result in an increased virial mass by a factor of ~2, 
furthermore optically thick, lower density J =1—0 transi¬ 
tions could cause an increase larger than a factor of ~5 if 
linewidths are more than double those measured here. It is 
therefore important that confirmed, optically thin tracers 
are used when calculating virial masses. 

The C^®0 gas masses in the sample are closely matched 
with the virial masses and are distributed about the 1:1 
line (Figure [TT] dot-dashed line), the results are consistent 
wit h those found for a smaller sample of IR-bright sources 
by iLopez-Sepulcre et al.l (l2010l i. However, the bisector fit 
to the data shows that the lower mass cores are actually 
skewed towards being unbound whereas the more massive 
ones are ~virialised, more clearly shown in Figure fT^ plot¬ 
ting the virtial ratio ve rsus core mass. This trend is seen by 
dUrauhart et al.ll2014al . and iKauffmann et al.l[201^ for indi¬ 
vidual datasets) where Myir/Mcore incre ases with decreas¬ 
ing clump mass. IKauffmann et aP (120131 ') argue that these 
trends can be explained if higher mass cores collapse and 
evolve rapidly through to the formation of their final stars, 
whereas lower mass cores may still have support present 
from, for example, outflow activity. Our results here are con¬ 
sistent with this picture. Only higher spatial resolution data 
will finally allow us to determine whether these global trends 
are reflected for individual protostellar sites, and how these 
interact with each other as a core collapses as a whole. 

There is no correlation between the surface density and 
radius, as expected for the Larson-like mass-radius scaling 
we observe in our data. There is however a scatter of almost 
an order of magnitude in the surface density. This correlates 
positively with the opacity in the line, in the sense 

that the highest surface densities have the highest opacities. 
However this can also be viewed in terms of the discussion re¬ 
garding whether the Larson relations are actually a function 
of the limiting surface density (or column density, opacity or 
extinction) seen for molecular cloud samples rat her than an 
underlying self-similar scal i ng relationship (e.g. Lada et "ni 
l201(j:lL ornbardi et aI.ll20 1Cll : iHeiderman et al.ll2010h . though 
see iBurkert &: Hartmannl (l2013l f for a discussion of pos¬ 
sible biases in the actual observational evidence for this. 
The rationale is that below this threshold there is insuffi¬ 
cient shielding to allo w molecular gas to form efficiently (see 
also the discussion in lEvans et al. I20l4l . The limiting star- 
forma tion surface density found by, e.g. IHeiderman et ahl 
(l201(tl . lies below the threshold of our sample as expected 
given Figure [HI We see no trend in the star-formation-rate 
surface density with gas surface density (as is evident from 
Mass-Luminosity plots, which essentially show the same 
data given the correct mass-radius correlation). However the 
correlation of surface density with opacity is consistent with 
this general concept of a threshold. 

iHever et al.1 (|2009l l pointed out that, if virial motions 
were the primary drivers of the Larson relationships, thus 
leading to a single surface density value regardless of size of 
cloud, then we should expect to see the ratio o-ul\fK as a 
constant. Figure [T^ shows the equivalent plot for our data. 
This clearly shows that this scaling depends on surface den- 
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Figure 11. Core mass versus the virial mass for all flag 0, 1 
(circles), 2 (squares) and 4 (triangles) cores. The open and filled 
symbols represent MYSOs and Hii regions, respectively. The plot¬ 
ted dot-dashed line is that of equal mass, while the dashed line 
indicates the bisector best fit. Most sources have Mcore ~ Mvir 
and are '^virialised, with a weak tendency for the lower mass 
cores to be less so than high mass cores. A 50 percent represen¬ 
tative uncertainty in core mass and an estimated uncertainty of 
15 percent in Virial mass (propagated from radius and FWHM 
uncertainties) are indicated to the bottom right. 



Figure 12. Virial ratio, Mvir/Mcore versus the core mass itself 
for flag 0, 1 (circles), 2 (squares) and 4 (triangles) cores. It is clear 
to see the more massive cores are closer to equilibrium, whereas 
the smaller ones have larger virial masses than the core mass. 
Representative 50 percent errors in the core mass and virial ratio 
(dominated by core mass uncertainty) are indicated to the bottom 
left. 
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Figure 13. The virial-surface density relation from iHever et all 
11^0 S '), for flag 0 and 1 (circles), 2 (sq uares) and 4 (trian - 
gles) cores. The crosses are taken from iHever et al.l |200j. 
The dotted line is that expected for virialised bound clouds, 
and the dashed line that for c louds in free-fall as outlined in 
iBallesteros-Paredcs et ^ |2oO). The propagated uncertainties 
are ~ 55 percent and ~ 7 percent for surface density and velocity 
dispersion/radius^/^ respectively. 


sity, just as they found (an d as seen in many other sample s 
including that of Larson — iBallesteros-Paredes et al.ll201l]) . 
In part this is consistent with the other results we have al¬ 
ready discussed. The scaling apparent in Larson’s relations 
depends critically on both the objects observed (i.e. real 
physical differences) and on the methods used. Crucially, the 
area over which the surface density is calculated should be 
matched with the appropriate kinematic data, and not cross 
matched with other measures, and as we have already noted, 
derived from an optically thin tracer that is representative of 
the structures in question. In principle if these simple guide¬ 
lines are followed we can compare relatively dis similar sam¬ 
ples. F igure 1 131 shows the s ample of clouds fromiHe ver et ahl 
ll2009l) . as well as our own. [Ballesteros-Paredes et al.l ( '2OIII ) 
also showed a s imilar plot for data from the infrared dark 
cloud sample of lCibson et al.l (l2009l ). 

We show bo t h the virial line expected for this plot 
from Hever et al. ||2009| ) and t he free-fall prediction from 
iBallesteros-Paredes et al.l ll201ll) . The data are best matched 
by the latter. This relationship is essentially carrying the 
same messages as the fact that mass, radius and FWHM 
form a fundamental plane, with the large velocity dispersion 
objects being those which are most massive for a particu¬ 
lar radius. It also essentially embeds the same result seen 
in Figure 1121 since lower virial ratios for more massive ob¬ 
jects tend to tally with higher velocity dispersions for higher 
surface densities as well. 

Finally we note that we find a strong linear correla¬ 
tion between the star formation rate su rface density (derived 
from the luminosity dKennicuttjl 19981 )) and the gas surface 
density ratioed with the free fall time. This is in agreement 


with the discussion in iKrumholz et ^ ll2012l ). but we would 
note that the correlation seen is inevitable given the form of 
the mass radius relationship. 


6 STAR FORMATION EFFICIENCY AND 
EVOLUTION 

6.1 Gas Mass and Luminosity 

Relationships between mass and luminosity can provide in¬ 
sights into star-formation efficiency and evolution in the 
cores. In our case the gas masses and source luminosities are 
independently established, whereas other works often obtain 
them via the same means, e.g. points in the mm/sub-mm 
SED. The downside to this is that many of the luminosities 
we have derived are at high spatial resolution, generally of 
the dominant source present within the original MSX beam. 
Therefore the C^®0 maps for such sources have effective 
radii larger than the “typical” beamsize for the SED from 
which we derive the luminosity. However, it is important to 
confirm that relationships found between continuum mass 
and luminosities are not due to effectively com paring the 
same data with itself (e.g. iMolinari et al.l(20081 ). Figure [TT] 
plots RMS source luminosity against core gas mass for both 
MYSOs and Hll regions (open and filled symbols, respec¬ 
tively). A linear relationship is observed for both MYSOs 
and Hll regions when plotting sources with mass flags 0, 1 
and 2 (and when including flag 4 sources). The fitted mass- 
luminosity relationship for MYSOs only is L oc 7V/i 08±0-50 
and for Hll regions only it is L oc ]\/[^- 0 ^roA 6 bisec¬ 

tor fit that best represents the scatter in the data for flag 
0, 1, 2 and 4 sources, a s the latter follow the same trend). 
lUrauhart et al.l ll2014bl ) find this slope continues for ~2 or¬ 
ders of magnitude in mass moving to much larger clump 
scales (Radius > Ipc). 


6.1.1 Probing core evolution 

Evolu tionary pre-main sequence tracks from Molinari et al.l 
(l2008l) , based on the model of iMcKee fc TanI 1 20031 ) , are 
also shown in Figure 1141 The basic physics that under¬ 
pins these tracks is relatively simple. Protostars continu¬ 
ously increase in luminosity, both through their strong ac¬ 
cretion phases and during Kelvin-Helmholtz (KH) contrac¬ 
tion, until they reach their even tual zero age main sequence 
(ZAMS) configuration (see, e.g. Hosolawa fc Omukail [20091 : 
iHosokawa et ai.ll2010l : IZhang et al.ll2014l ). Once on the main 
sequence, massive and intermediate mass stars fairly quickly 
disperse their natal molecular material through the action 
of their ionising radiation and winds. 

We have also plotted the location of the stellar 
ZAMS in Figur e 1141 We take stellar luminosities from 
ISalaris fc Cassisil (l2006l . Figure 5 .11) for mas s es ran ^in^ 
from 0.5 Mq to 6 Mq and from IPavies et al.l (l201ll ) for 
masses >6M0. We assume that the total mass of stars 
Afstars = Mcoie and the protostellar masses are distributed 
according to the Salpeter power-law IMF (using only stellar 
masses from 0.5 to 150 Mq). This assumption is equivalent 
to a star-formation effi ciency (SFE) of 50 percent, follow- 
ing lLada fc Ladal ||2003| ) [SFE = Mstars/(A/stars + ALgas)], if 
no gas has been lost due to winds or outflows. This ZAMS 
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Core Mass 


Figure 14. Mass-luminosity plot for all MYSOs and Hll regions, with mass flags 0 and 1 as circles, flag 2 sources as squares and 
flag 4 sources as triangles. Open and filled symbols separate MYSOs and Hn regions. The two star symbols are the flag 7 Hll regions 
G077.9637—00.0075 and G080.9383—00.1268 that appear to have dispersed their core material (assuming the cores were morphologically 
similar to the other cores in the sample at an earlier evolutionary stage). The dot-dashed line is the ZAMS luminosity expected from the 
most massive star in the core, when the stellar masses are distributed according to the Salpeter IMF in the case of a 50 percent SFE 
(the mass in stars equals the core mass). For the ZAMS line to fit the data, higher SFEs are required (moving it upwards). The dashed 
line is the luminosity of a cluster of ZAMS stars distributed with a Salpeter IMF with a 50 percent SFE. A reduced SFE is required if 
one were to use the luminosity from the whole cluster as it is already overestimates the lumi nosity of the m o st ma ssive cores (see Section 
16.1.211 . The dotted lines represent the pre-main sequence tracks extracted from models bv iMolinari et all ll2008h and can be used as a 
tool to trace evolution. Illustrative error bars of 50 percent in mass and 30 percent in luminosity are in the right-hand bottom corner. 


line is shown as the dot-dashed line on Figure 1141 and the 
evolutionary tracks as dotted lines, for initial core masses 
of 80, 140, 350, 700 and 2000 Mq . The turnover in the 
IMolinari et al.l ll2008h tracks towards envelope dispersal oc¬ 
curs on this ZAMS line as expected. A change in the SFE 
is equivalent to shifting this ZAMS line right or left (the 
luminosity does not change, but the residual core mass is 
either greater for a lower SFE or less for a greater SFE). 
This is discussed further in the next section. We have not 
attempted to model a cluster of stars in detail, given that 
high mass stars evolve more quickly than lower mass stars, 
and hence reach their ZAMS luminosity more quickly. At 
least for clusters of protostars this should mean that the 
higher mass members do minate the luminosity (the same is 
not true for Hll regions - iLumsden et al.ll2003l) . 

The relatively small scatter in this diagram suggests 
that most of the objects we have observed are of s imilar 
evolutionary stage, as noted bv lUrauhart et ^ ll2014bl) . It is 
worth noting that the Hll regions overlap with the MYSOs, 
which suggests that our ‘compact’ size criterion ensures 


that for th e most part we have selected relatively young 
Hii regions. IMolinari et al.l (l2008l ) note that the model PMS 
tracks are variable and dependent upon accretion and core¬ 
dispersal rates, for example, such that they can shift both 
vertically and horizontally which may partially explain this 
spread, as could the observational errors in both mass and 
luminosity. Furthermore, these tracks assume that the core 
does not gain mass from its surroundings by large-scale in¬ 
fall during star formation, and that one core leads to one 
star with a fixed star formation efficiency. If, for example, 
a core gains mass during the formation process, the tracks 
would also slant to the right as they increase in luminosity 
as they gain extra mass. An alternative explanation however 
is given by the timescales involved for each phase of evolu¬ 
tion as a function of final mass. The most massive stars reach 
the ZAMS (and hence power Hll reg ions) much more rapidl y 
than it takes a less massive star (e.g lUrguhart et al.ll2014^ . 
Therefore the most massive stars reach the ZAMS and form 
Hii regions whilst still very heavily embedded. Notably, the 
Hii regions of lower mass appear more evolved on average 
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in this figure than those of higher mass, which suggests that 
this is what we are seeing. Although it is likely all scenarios 
contribute to some extent. The filled star symbols in Figure 
1 141 show the two Hll regions flagged as 7 that are thought to 
be the most evolved and have dispersed their cores. Their 
position to the far left hand side of the mass-luminosity di¬ 
agram su pports this inter p retatio n. _ _ 

Both iMottram et al.l (l2011ah and iDavies et al.l ll201lh 
find the lifetimes of MYSO and Hii-region phases to be 
roughly comparable fof the order 10® yr). In particular, Fig¬ 
ure 7 of I Davies et al.l (l201ll ') illustrates how the source clas¬ 
sifications for stars with final masses >8Mq change (from 
MYSO to Hll region) over a very narrow period in time 
where core masses and luminosities are still coincident. As 
this paper is a pre-cursor to a detailed investigation of molec¬ 
ular bipolar outflows, a narrow distribution of evolutionary 
stage makes this sample ideal for identifying trends in out¬ 
flow parameters due to source properties, free from any ef- 

fect s of a spread in s ourc e ages. _ 

lElia et ^ (l2010l l and iBeltran et all il2013l l use the L/M 
ratio as an evolutionary tracer, where more evolved sources 
are more luminous and have dispersed more material from 
their cores/envelopes (increasing L/M) in comparison with 
lower luminosity, deeply embedd ed sources (th e ‘evol ved’ 
flag 7 sources have L/M > 10®). iBeltran et al.l (l2013tl in¬ 
dicate a clear dichotomy between 24 ^m-bright and dark 
sources in their survey of the G29.96-0.02 cloud. In a sim¬ 
ilar vein, the distribution of L/M in the Hll regions are 
examined, looking for differences with MYSOs. Of course 
one caveat is that we assume that to produce a source of 
a certain luminosity we always start at the same mass (i.e. 
star formation efficiency does not vary). It is not clear that 
this is the case from core to core, and therefore can introduce 
some scatter in the distribution of sources (as previously dis¬ 
cussed). Figure fTSl shows the L/M histogram for MYSOs and 
Hll regions, and shows no clear difference between the two 
subsamples. Both a Kolmogorov-Smirnov (KS) test and a 
Mann-Whitney U (MWU) test indicate that MYSOs and Hii 
re gions are drawn from th e same distribution (as also found 
by Urauhart et al.|[2014bll . The reported D value (~0.17) is 
not above that required (>0.32) for the sources to be de- 
tectably drawn from different distributions in the KS test. 
The reported MWU test probability (~0.46) is well above 
the 0.05 significance level and so the hypothesis that the 
two samples are drawn from the same distribution is not 
rejected. 


6.1.2 Star formation efficiency from cores to stars 

The masses plotted in Figure [14] are those of the clumps 
containing what is probably a cluster of protostars within 
which the massive protostar detected by the RMS survey 
is forming. However, the luminosities could either be dom¬ 
inated by the most massive star, or be the total of all the 
protostars in the cluster. The dot-dashed line in Figure [14] 
represents the luminosity of a single star, the most massive 
in the cluster. In principle, this line should set an upper en¬ 
velope to the observed data for single stars, and allow us 
to derive an estimate of the star formation efficiency (SFE). 
As noted above, although clusters of protostars are probably 
dominated by their most massive member, the cluster line 
still lies somewhat higher in luminosity than the single star 


10 


8 



0.5 1.0 1.5 2.0 2.5 

Luminosity/Mass LogfL^/Ms) 


Figure 15. Histogram of the Luminosity/Mass ratio for MYSOs 
and Hii regions, respectively, in light and dark grey bars. Both 
a KS and Mann-Whitney U test cannot distinguish MYSOs and 
Hii regions as different distributions. 


line (dashed line). The shape of the single star line matches 
the data reasonably well, suggesting that this simple model 
is reasonable, i.e. the most massive star dominates. The clus¬ 
ter line however, is steeper. This is likely due to the fact that 
we plot ZAMS luminosities, which is a reasonable assump¬ 
tion for the more massive protostars, but probably underes¬ 
timates the luminosity of accreting low mass protostars in 
theses clusters. 

The overall SFE should be the value that sets all of 
our observations below the actual model lines. If each of our 
objects genuinely only contained a single luminous source, 
or else they are on the track beyond the ZAMS and into 
the core dispersal region, they will fall below the dot-dashed 
line. This requires a shift to lower core masses (higher effi¬ 
ciencies) by close to a factor of two. Such a large shift, and 
high SFE (>60percent), appears implausible. Instead the 
alternative that we are seeing clusters (dashed line), still in 
agreement with our observations, which enhances the total 
luminosity at larger core masses by about a factor of two, 
seems plausible. In this case the SFE would be in the 40- 
50 percent range. 


7 SUMMARY 

The core parameters of 94 sources (of 99) MYSOs and Hii 
regions selected as outflow candidates representative of the 
pre-2008 RMS survey, have been established. The sample is 
reduced to 89 when enforcing the distance limit of 6 kpc and 
is still representative of massive protostars across the MYSO 
to Hii region transition in the up-to-date RMS survey that 
now meet the original selection criteria. 

The majority of the cores exhibit a single, Gaus¬ 
sian G'^^O line profile. Larson-like relationships (FWHM 
linewidth relationships with luminosity, radius and mass) 
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are found for all the cores, though the scaling for these is 
not continuous with other studies of larger-scale molecu¬ 
lar clouds. The most fundamental relationship is between 
mass and core radius, which gives rise to a surface den¬ 
sity independent of radius. We note, however, that the scat¬ 
ter in the surface density in this relation is well explained 
by the correlation with gas opacity and agrees with mod¬ 
els in which Larson-style relations arise due to the ob¬ 
served surface or column density limits. We find two possi¬ 
ble fundamental planes in this work, representing core evo¬ 
lution (mass-luminosity-radius) and likely virial contribu¬ 
tions (mass-radius-FWHM). All cores appear to be virialised 
Mcore ~ Afvir. Core parameters on the observed scales are 
interpreted as being set prior to the onset of star formation 
and are not subsequently effected by the feedback from the 
massive protostars. 

A correlation is found between the dust-continuum 
masses and the core (gas) masses established from the C^®0 
emission. At the resolution of the observations, both the 
dust-continuum emission and C^®0 emission trace the same 
material and structures in both MYSO and Hll regions. The 
small differences found are all consistent with dust being 
a better tracer of diffuse low density gas on larger scales, 
whereas the C^®0 traces the dense hearts of the molecular 
clumps. 

The sources are consistent with most of the luminosity 
arising in a single massive protostar (or central ZAMS stars 
for the Hll regions). The tight banding of all sources in the 
M-L plot, accounting for scatter, indicates a similar evolu¬ 
tionary stage for both the MYSOs and Hll regions investi¬ 
gated. This means they are ideal candidates to investigate 
relationships between outflow and source properties between 
the source types to be detailed in an upcoming paper. 

Further development of the RMS survey database is 
required with higher spatial resolution at millimetre/sub¬ 
millimetre wavelengths to understand how the core masses 
are specifically related to the immediate regions surround¬ 
ing the most massive protostars on arcsecond scales. This 
will allow us to investigate whether the mass distributions 
are different for MYSOs and Hll regions, if CO depletion is 
present in these cores and on what spatial scale, whether 
cores have further velocity sub-structure influenced by out¬ 
flows, and how Larson-type relationships hold below 0.1 pc 
scales. 
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APPENDIX A: TEMPERATURE, OPTICAL DEPTH, COLUMN DENSITY AND MASS 
CALCULATIONS 

In this appendix the column density and mass equations are derived following from the result of iGarden et alj (llQQlTl . except 
for the CO(3—2) transition. The total column density of a linear, rigid rotor molecule under conditions of local thermodynamic 
equilibrium (LTE), with the populations of all levels characterised by a single excitation temperature, Tex, is obtained from 
the integral of the optical depth over the line profile: 

_ 3fc exp[/ii?J( J + l)/ferex] Tex + hB/3k 
8 'K^Bfi'^ (J +1) [1 — &q){—hv / kTe:^)] 

where B is the rotational constant, /r is the permanent dipole moment of the molecule and J is the rotational quantum 
number of the lower state, in this case J = 2 for the CO(3—2) transition, k and h are the Boltzmann and Planck constants 
respectively. 

Both the excitation temperature Tex and optical depth ris are solved for. Assuming both the source function and initial 
intensity are blackbodies at the respective temperatures Tex and Tcmb = 2.73 K, then 


J- 


, dv 


(Al) 


Tmb = 


hv 


1 

exp(/iiz/A:rex) — 1 


1 

exp{hiy/kTcmh) - 1 


X [l-exp(-T)] , 


(A2) 


for any line transition. From Equation I A2I the excitation temperature, Tex, can be derived, provided the optical depth of the 
line is known. The *^^00 observed simultaneously with C^®0 should be optically thick at the locations of peak emission in 
our cores, such that as T13 00, [1 — exp(—r)] —>■ 1. In the case that 0*^*0 is optically thin, we first calculate the optical 

depth of ri3 directly from the ratio of the observed and C^®0 antenna temperatures, in order to confirm is thick: 


r,nb,i3 _ 1 - ^ 1 - 

Tmb.is “ 1 - e-^18 1 - e(-^i3/R) 

where R is abundance ratio of [^^CO]/[C^®0] derived from. 


(A3) 


[^^0]/P^0] ^ 58.8 xD,,+37.1 

[i2C]/[i3C] 7.5xDgc + 7.6 '' '' 

following IWilson fc RoodI lll994l and Dgc is Galactocentric distance in kpc. 

The excitation temperature of *^^00 is now calculated. Rearranging Equation IA2I and substituting in the numerically 
established optical depth, ri3. 


Te, 


15.86 

In [1 + 15.86/ ({r,nb,13/(1 - exp(-ri3))} + 0.047)] ’ 


(A5) 


where hv{}^CO)/k = 15.86 K, with (/('^^GO) = 330.58 GHz and Tmb.is is the main-beam brightness temperature of the ^®GO 
emission. 

The 0*^*0 excitation temperature Tex is assumed to be the same as that of the line, and thus the optical depth, 

ri8 can now be calculated for the 0*^*0 emission from: 


ri8 = -In 


1 - 


Tmb.l 


15.80/[exp(15.80/Tex) - 1] - 0.045_ 

Here the approximation for J dv follows feuckle et al.l ll2010h for the case where r y^O: 


(A6) 


/ 


dv = 


hv f 1 

k \exp(/iiz/fcTex) — 1 


_ 1 _ 

exp(hi//fcTcmb) - 1 


-1 


T 

[1 - exp(-r)] 


/ 


Tmb dv 


(A7) 


The brightness temperature, Tmb, is the antenna temperature of the telescope divided by the beam efficiency, TX/iimb, 
and corresponds to the Rayleigh-Jeans brightness of a source minus the brightness of the cosmic microwave background with 
temperature, Tcmb = 2.73 K, over the beam. Combining equations ED and EH in the limit where Tex ^ Tcmb results in the 
column density: 


3k exp[/iBJ(J-I-l)/fcTex] 1 Te^ + hB/3k 7 1 , . 

8 -k^Bij? (J-fl) (hiz/fe) [exp(-hi//fcTex)] y [1 - exp(-r)] ^ ’ 

where the permanent dipole moment for is 0.1101 Debye (IChackerian fc TipDind[i983fl . Conforming to cgs units typically 
used in such analysis i?=58.14GHz, fc=1.381xlO”*^® ergK“^, h=6.626x 10~^^ ergs, iz(C'^®0) = 329.33055 GHz, velocity v is in 
kms“^, ^(C^®0) = 0.1101 X lO”*^® StatC cm (where 1 statC = 1 g^^^ cvo?^^ s^ = 1 erg^^^ cm^^^), r becomes ns, the calculated 
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Figure Al. Correction factor between calculated LTE calculated column density and that input into RADEX. Although at 20 K the 
correction factor is ~2 at 10'*cm“^, at the average volume density of the cores (~7xl0‘^ cm“®) the correction factor is ~ unity even at 


20 K. 


optical depth of the C^®0 line calcualted in Eouation lAdl and Tex is the calculated excitation temperature from Equation IA5I 
The column density for the C^®0 (3—2) transition is therefore: 


Af(C^®0) =5.0x10^^ 


exp(16.74/rex) (Tex + 0.93) 
exp(-15.80/Tex) 



r-is 

[1 - exp(-Ti8)] 


dv , cm 


The mass of a source can then be calculated from the column density via: 


(A9) 


ng D'^ (AlO) 

where fig = 1.36 is the total gas mass relative to H2, the ab undance ratio [H2/C^®0] is a combination of H2/^^CO = 10^ and 
160/I80 = 58.8 X DGc(kpc) + 37.1 l| Wilson fc Roodlll994l l. where Dqc is the Galactocentric distance and D is the distance 
of the source to the Sun, both in kpc. 11 is the solid angle corresponding to the emission in one pixel of the maps used in this 
work. Thus including the conversion factors the core gas mass in solar masses (M0) is calculated for every pixel of the map 
using Equation lAl II The total core mass is the summation within the defined aperture encompassing the core, as described 
in Section [3] 


Mg,, = Af(CO) 


H2 


CiSQ 


Mg„ (Mq) = 2.5 X 10“^^ e^C) D\kpc) 


H2 

CiSQ 


exp(16.74/rex) (Tex + 0.93) 
exp(-15.80/rex) 



ns _ 

[1 - exp(-ri8)] 


dv 


(All) 


The LTE approximation is only valid if the volume density is sufficient enough (n>nci.it = ~3.2x lO'^ cm“®) 0 such that the 
excitation temperature can be assumed to be equal to the kinetic temperature. If this is not the case (Tex < Tkin) the LTE 
calculation will underestimate the column density and therefore mass. These parameters are then subject to a correction 
factor dependent on excitation temperature and volume density. Figure shows the correction factor as a function of 
volume density for kinetic temperatures of 20, 30 and 40K (i.e. the input excit ation temperatures in LT E calculations). The 
factor is calculated by comparison of the input column density into radex (Ivan der Tak et ^1200711 and that calculated 
via LTE calculation here. A representative H2 column density of 2xl0^^cm“^ and a linewidth of the average of all cores 
are used as inputs. At 20 K for a 10"^ cm“® volume density (just below critical), the correction factor required to obtain the 
correct column density from the presented LTE assumption is ~2. Above the critical density the correction factor approaches 
unity. The average volume density of the cores however is ~7xl0'’cm“® calculated using the LTE masses and deconvolved 
radii. Thus the LTE assumption is valid as the correction factor at this volume density is ~ unity. There is no evidence for 
sub-thermal excitation of the CO J=3—2 transition. 


® see, http://home.strw.leidenuniv.nl/~moldata/ 
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